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Recently we proposed a particle-number-conserving theory for nuclear pairing [Jia, Phys. Rev. 

C 88, 044303 (2013)] through the generalized density matrix formalism. The relevant equations 
were solved for the case when each single-particle level has a distinct set of quantum numbers and 
could only pair with its time-reversed partner (BCS-type Hamiltonian). In this work we consider 
the more general situation when several single-particle levels could have the same set of quantum 
numbers and pairing among these levels is allowed (HFB-type Hamiltonian). The pair condensate 
wavefunction (the HFB wavefunction projected onto good particle number) is determined by the 
equations of motion for density matrix operators instead of the variation principle. The theory is 
tested in the simple two-level model with factorizable pairing interactions and the semi-realistic 
model with the zero-range delta interaction. 
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I. INTRODUCTION 

The Bardeen-Cooper-SchriefFer (BCS) theory [l| or 
its advanced version the Hartree-Fock-Bogoliubov (HFB) 
theory 0 has long been used to treat pairing correlations 
Q in atomic nuclei. In the mean field we introduce the 
Bogoliubov quasi-particles and write the ground state as 
a Slater determinant of the latter. Usually the variation 
principle is used to determine the structure of the quasi¬ 
particles. Although enjoying great success, the method 
has disadvantages of breaking the exact particle num¬ 
ber and a need for an unphysical minimum pairing-force 
strength Si. 

The theory could be improved by the “variation af¬ 
ter particle-number projection” (VAPNP) procedure, in 
which the BCS or HFB wavefunction is projected onto 
good particle number before the variation is done @ . Ef¬ 
fectively, the pair condensate wavefunction [Eq. Q] is 
taken as the variational ground state. The method has 
been discussed [Ml and applied with success to ultra¬ 
small metallic grains |lll - H3 i n the VAPNP+BCS ver¬ 
sion, and to realistic nuclei [l4l - Ha| in the VAPNP-bHFB 
version. However, the computing time cost is relatively 
large and there is only a limited number of nuclei on the 
nuclear chart that has been calculated by the VAPNP -b 
HFB method in large configuration spaces. 

Recently we proposed a new criteria to determine 
the pair condensate wavefunction based on the Heisen¬ 
berg equations of motion for density matrix operators. 
The relevant equations have been solved for the BCS- 
type Hamiltonian: Each single-particle level has a dis¬ 
tinct set of quantum numbers corresponding to the sym¬ 
metries of the pairing operator, thus a level could only 
be paired with its time-reversed partner. The validity 
of the theory has been proved on a large ensemble with 
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random interactions [2^. In this work we consider the 
general situation of the HFB-type Hamiltonian, where 
several single-particle levels could have the same set of 
quantum numbers and pairing among them is allowed. 
The method should cost much less time in computing 
compared with the traditional variational calculation. In 
the algorithm of the current theory only the one-body 
density matrices on the pair condensate are calculated, 
whereas the variation principle needs the two-body den¬ 
sity matrices in computing the expectation value of a 
two-body Hamiltonian. 

In Sec. miwe give the general formalism to solve for the 
particle-number-projected HFB wavefunction with the 
pairing theory based on the generalized density matrix 
(GDM) formalism. Then the method is tested in the 
simple two-level model with factorizable pairing interac¬ 
tions and in the semi-realistic model with the zero-range 
delta interaction, in Sec. m and Sec. IlYl respectively. 
Finally Sec. IVl summarizes the work. 

II. FORMALISM 

The antisymmetrized (two-body) fermionic Hamilto¬ 
nian governing the dynamics of the system is written as 

H — ^ ^ -b “ ^ ^ (1) 

Oi0 a(3-y5 

In the presence of pairing correlations, the ground state 
of the 2iV-particle system is assumed to be an V-pair 
condensate, 

= ^iPY\o), ( 2 ) 

yJXN 

where xn is the normalization factor, is the pair cre¬ 
ation operator 

aGA,l3^A 
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In Eq. dS]) A represents a subspace consisting of “pair- 
indices” whose dimension is half of that of the single¬ 
particle space. We could take, for example, A to consist 
of those single-particle levels with a positive magnetic 
quantum number. |/3) is the time-reversed level of the 

single-particle level |/3) {[P) = —1/3)). The pair struc¬ 
ture matrix Vap is Hermitian and block-diagonal. The 
time-reversal invariance of \(j)N) implies that v^a = 

(Tua/ja^ at|0) = Vap vanishes unless \a) 

and 1/3) have the same set r of quantum numbers related 
to the symmetries of Pb For example, in the spherical 
shell model t = (tt, j, m) is a collection of parity tt, angu¬ 
lar momentum j and its projection m; in the deformed 
Nilsson mean field r = (tt, K) is a collection of parity tt 
and angular-momentum projection K onto the symmetry 
axis. 

We introduce the unitary transformation ry between 
the original single-particle basis {la)} and the new single¬ 
particle basis {|1)} as 

Vai = (a|l), (4) 

where elements of the transformation matrix ry are de¬ 
fined through Dirac notation. Properties of time-reversal 
operation imply that (ry^jia = (l|a) = (a|l). The matrix 
7y is block-diagonal: Tyai vanishes unless |1) and |a) have 
the same quantum number r. Under the transformation 
(HI) the operator ([31) becomes 

pt = ^ ui2a|oi, (5) 

leA,2GA 

with 

Vl2 = ^ {ri^)laVaPVl32- (6) 

aGA,/3eA 

We choose the transformation ry that diagonalizes the 
Hermitian matrix Vap in Eq. ®, consequently be¬ 
comes 

pt=^^,la^a^. (7) 

leA 

We call this new single-particle basis the “canonical ba¬ 
sis”. In the following the Arabic numerals |1), |2), ... 
refer to single-particle levels in this basis unless other¬ 
wise specified. 

In the canonical basis the density matrices for the pair 
condensate m are “diagonal”: 

Pi2 = {(j)N\alai\(j)N) = Si2ni, (8) 

ki 2 = (4>N-i\a2ai\4'N) = S 12 S 1 . ( 9 ) 

The normalization factor xn occupation numbers ni 
([8|). and pair-transfer amplitudes si m are functions of 
the pair structures vi ©; their functional forms, as the 
“kinematics” of the system, have already been given in 
Eqs. (23) and (24) of Ref. [T^ and are not repeated here. 


The Hartree-Fock mean held / and pairing mean held 5 
are dehned as 

/12 = ei2 + R1432P34 = ei2 + Ri332’T-3j (10) 

34 3 

'^12 = 2 ^ ^1234^43 = 2 ^ lb233®3- (H) 

34 3 

/ and 5 are block-diagonal matrix: /12 and 5^2 vanish 
unless |1) and |2) have the same quantum number r. 
The Hamiltonian parameters (ei2 and V1234) in Eqs. (flUl) 
and (EH) should be calculated from those in the original 
single-particle basis {tap and Vap^ys) through the trans¬ 
formation Q. 

The equation of motion for the density matrix k has 
been derived in Eq. (14) of Ref. (iH| . 

{En - E]sr^i)K = fK +Kf'^ + S - Sp^ - pS, (12) 

where E^ and Pat-i are ground state energies for \(j)N) 
and \(j)N-i), respectively. On the right-hand side, / and 
6 are the mean helds Eon and EH, is the transpose 
of /, and terms like ‘/k’ are understood as matrix mul¬ 
tiplication. In deriving Eq. EH we have used the main 
approximation of the method, 

{(j>N-i\a\a3a2ai\(j>N) 

Ri {(j)N\alai\4>N){4>N-i\a3a2\(l>N) 
-{(j)N\ala2\4’N){4’N-i\a3ai\(l>N) 
+{(j)N\ala3\(j)N){4’N-i\a2ai\(l>N) 

= P 14 K 23 — P 24 K 13 -I-P 34 K 12 , (13) 

which says that on the pair condensate 0 the two- 
body density matrix factorizes into products of one-body 
density matrices in both the particle-hole and particle- 
particle channels. 

Equation (fl^ is a block-diagonal matrix equation; its 
12 matrix element vanishes unless |1) and |2) have the 
same quantum number r. Within each block we take the 
12 matrix element on both sides of Eq. EH, when 1 2 

we get after simplification 

0 = (si-I-S2)/i2 + (1 - - n2)(5i2- (14) 

And the 11 matrix element gives 

En — En-i = 2fii + 6^1 -( 15 ) 

Si 

Equations EH and (fT51) are the main equations of the 
theory. Below we show that the number of constraints 
from these two equations equals to the number of param¬ 
eters in the pair structure Vap 0; thus the latter is fixed 
completely. 

We assume a single-particle space of dimension 2D split 
by the quantum number r as D = Qt- The number of 
restrictions in Eq. (HH is Er = E,- llr (fir - 1) = 
Et — D (there is a factor 2 because the equation has 
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real and imaginary parts). Equation (jl5l) implies that 
the right-hand side is independent of the single-particle 
label 1, which gives — 1 constraints. Hence the total 
number of constraints is ~ other hand, 

the number of independent parameters in the Hermitian 
block-diagonal pair-structure matrix Vap (tS]) is ~ 1 

(There is a “—1” because an overall normalization fac¬ 
tor does not matter). The number of restrictions indeed 
equals to the number of parameters. 

In practical calculations usually the Hamiltonian pa¬ 
rameters eap and Vap^s o are real. In this case the 
pair structures Vap could be taken as real numbers 
and the transformation rj @ is an orthogonal matrix. 
The number of restrictions from Eqs. o and m is 
-I- H — 1 = — 1, which equals to the 

number of independent parameters in the real symmet¬ 
ric block-diagonal pair-structure matrix Vap ([3]). In the 
following we assume that this is the case. 

An alternative parametrization may be more conve¬ 
nient in practice. The independent parameters could be 
taken as those in the orthogonal transformation ry o 
and the pair structure in the canonical basis vi ([7]). The 
orthogonal transformation for H,- = 2 could be parame¬ 
terized as 


V = 


For = 3, 



ria2 ' 

1 ( COS 0 

— sin 0 \ 

V ^/31 

VP2 ) 

' ~ sin0 

COS0 y 


(16) 


/ cos 013 

0 

-sin 013 \ 

, /I 

0 

0 \ 

0 

1 

0 

0 

cos 023 

-sin 023 

V sin 013 

0 

cos 013 j 

1 0 

sin 023 

cos 023 / 


And in general, 


7] = _[J_77 
i<j 




(17) 


III. TWO-LEVEL MODEL 


We test the formalism in the simple model with two 
single-particle levels and a factorizable pairing interac¬ 
tion. The theory is solved analytically and the results 
are compared with the exact shell-model results. 

We assume the rotational invariance and a single¬ 
particle space of two j-levels each with degeneracy 20 = 
2j -\-1 (degenerate in the magnetic quantum number m). 
The fermionic Hamiltonian takes the form 




Oi^ Cirri'^ 


+ ^pO^pm^Pm) - n'l'H, (18) 




cos 012 —sin 012 0 
sin 012 cos 012 0 

0 0 1 


where H^ is the pairing operator 


^ ^ ^ id^am^ani A 9^pm^Prh P^am^prh P^Pm^aih)^ 

m>0 

in which g and p are “diagonal” and “off-diagonal” pairing strengths, respectively. The Hamiltonian parameters in 
the canonical single-particle basis are calculated through the transformation (IT6t . 

= ( 20 ) 

m ij 


where 


and 


( £11 £12 \ / £„ cos^ 0 -I- e/3 sin^ 0 i(e/3 - Cq.) sin 20 \ 

E21 ^22/ \ i(e/3 - Cq.) sin 20 Cq, sin^ 0 -I- e/3 cos^ 0 E 


( 21 ) 


= E [(5 + ^ sin 20)aEaE + (fl - P sin 20)aEaL + P cos 200^4™ + P cos 20aEah] ■ (22) 

m>0 

Based on Eqs. (HI and (HU we compute the matrix elements of the mean fields in the canonical single-particle basis, 

fii = £a cos^ 0 + ep sin^ 9 — cos^ 20n2 — {g + P sin 20)^ni, 

/22 = £a sin^ 0 -I- e^ cos^ & — [g — psin20)^n2 — cos^ 20ni, 

/12 = ^(£/3 ~ ea)sin20 — pcos20[(p — psin20)n2 -I- {g -|-psin20)ni], 
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and 

= —^{g + psin29)[{g +psin20)si + {g — psin20)s2], 

S 22 = —^{g sin 2d) [(5 +psin2d)si + {g — psin2d)s2], 

6^2 = —^pcos29[{g +psin2d)si + {g — psin2d)s2]- 

Consequently Eqs. (HI and (ITSl) imply 

0 = sin2d(si + S2)(e/3 - £„) + 2gpcos29{si + S2)[(»^i + n 2 ){^ - 1) - O] 
sin4d[12(ni + 712 — l)(si — S2) — (si + S2)(?^i — ^12)], 


(23) 


and 

0 = cos2dsiS2(ea — £/?) cos^ 2dsiS2(7ii — 712) + sin^ 2d)[(si)^(l — 2712) — (^2)^(1 — 27ii)] 

—2gp^su\29siS2 + (12 — l){g + psin2d)^7iiSiS2 — (12 — 1)((7 — psin2d)^r72SiS2- (24) 

These two equations fix the two independent model parameters V 2 /V 1 and 9 (ni and si are functions of V 2 /V 1 ). Then 
the density matrices in the original single-particle basis are calculated through the transformation g (11611 . 

Paa Pap\ _ ( 7li COS^ d -f 772 sin^ d (tIi — 7I2) Sin d COS d 
Ppa d/3/3 ) \ (tii — 772) sin 9 COS 9 m sin^ d -|- 772 cos^ 9 



and 


/ Kaa _ ( Si cos^ d-I-S2 sin^ d (si — S2) sindcosd 

\ K^a Kpp ) \ (si — S2) sindcosd Si sin^ d -I- S2 cos^ d 


(26) 


We consider the range of model parameters in the re¬ 
alistic spherical nuclear shell model. Usually the two 
single-particle levels ja) and |/3) belong to different major 
shells and ep — ta is bigger than 10 MeV in magnitude 
that is about the energy of two major-shell gaps (adja¬ 
cent major shells have opposite parity). Based on the 
empirical pairing strength formula G = g'^ Ki 20/A MeV 
{A is the mass number) [2l| , G = g'^ should be less than 2 
MeV in medium and heavy nuclei having considered pos¬ 
sible deviations from the constant pairing formula. The 
off-diagonal pairing strength p is usually smaller than the 
diagonal pairing strength g owing to the smaller overlap 
of the two single-particle wavefunctions involved. 

We test the model numerically in an ensemble consist¬ 
ing of examples with different model parameters. The 
single-particle angular momentum j of the examples 
takes value from j = | > f, || ■ The number of pairs 
is in the range 1 < iV < 2j = 2f2 — 1 (from two parti¬ 
cles to two holes in the model space). The energy gap 
between the two single-particle levels is selected to be 
the energy unit so that = —0.5 and ep = 0.5. Based 
on the estimations in the previous paragraph, we choose 
the range of the pairing strength to be 0 < g < 0.5 and 
0 < p < 5 with step size 0.1. We scan the whole pa¬ 
rameter space and the total number of examples in the 
ensemble is 15 2j = 15 x 24 = 360 (the factor of 15 is 
the number of possible values of g and p). 

The results for the density matrices p and k in the 


original single-particle basis are shown in Figs. [T] and [U 
respectively. Each point in these figures corresponds to 
one example in the ensemble with the horizontal coor¬ 
dinate being the exact result of the respective quantity 
and the vertical one being that from the GDM calcula¬ 
tion, thus a perfect calculation would have all the points 
lying on the y = x straight line. From Figs. [Hand[2] we 
see that in general the GDM theory reproduces the exact 
density matrices well. The root-mean-square deviations 
are 

O'Pec = = 0.0125, = cTpp^ = 0.0198, 

= 0.0183, = 0.0214, = 0.0384. 

o'pee and cFpfif, are equal because paa -f ppp = N/fl. 
The results for the ground state energies are shown in 
Fig. [31 The horizontal coordinate of each point is 
the pairing correlation energy for the example Fpair = 
£i 77 f — Fexact, where Fexact is the exact ground state 
energy by the direct diagonalization, and nf = 1 or 0 
is the occupation number of the naive Fermi distribu¬ 
tion. The vertical coordinate is the ground state energy 
by the GDM calculation measured from the exact one, 
Egum = {(I)n\H\(I)n) - Fexact, where is cal¬ 

culated in the canonical single-particle basis using the 
Hamiltonian dlO]) tog ether with the recursive formulas 
derived in Refs. 113, 123 ■ We see that in general the 
GDM calculation reproduces well the ground state en¬ 
ergies: the errors are small compared with the pairing 
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correlation energies. The average values of the errors 
and the pairing correlation energies for the ensemble are 
.Egdm =0.0198 and impair = 1-78, respectively. 

At last in Fig. 0]we show the accuracy of the method 
depending on the off-diagonal pairing strength p. The 
360 examples in the ensemble are divided into five sub¬ 
groups according to the value of p, and deviations within 
each subgroup of various quantities are calculated. We 
see that in general the errors of the GDM calculation in¬ 
crease with p. In realistic spherical nuclear shell model 
p should be smaller than about 0.3 based on the esti¬ 
mations made in the paragraph below Eq. (1261) . In this 
region the GDM theory is rather accurate. 

IV. MODEL WITH ZERO-RANGE DELTA 
INTERACTION 

In this section we further test the theory in a semi- 
realistic model with the pairing Hamiltonian matrix el¬ 
ements calculated from the zero-range delta interaction. 
The model has only one species of nucleons (for example, 
neutrons) and the single-particle space consists of five lev¬ 
els: three of them { 2 si/ 2 , 1 ^ 3 / 2 , l'^ 5 / 2 ) from the AHlu 
major shell and two (3si/2, 2 ^ 5 / 2 ) are from the 6huj major 
shell. The single-particle energies are determined by ran¬ 
dom number generator to be £ 281/2 = 1-16, £ 1 ^ 3/2 = 2.II, 
£id 5/2 = O.II, £ 3 si /2 = 12.98, and £ 2 ^ 5/2 = 10.37 MeV 
(The gap between the ifnj and Ghoj major shells is taken 
to be around 10 MeV). The two-body pairing matrix el¬ 
ements [(a/3; J\V\'-fS-, J) with J = 0] are calculated from 
the zero-range delta interaction V’(ri, r 2 ) = — Ai5(ri —r 2 ), 
taking the single-particle wavefunctions (|a), |/3) ...) to 
be the harmonic oscillator ones. We perform six sets of 
calculations at different model parameters of the parti¬ 
cle number 2N and the pairing strength A as listed in 
table m The unit for A is MeV • 5^, where b = ^/K/Jpio) 
is the length parameter for the harmonic oscillator wave- 
functions {pL is the mass of the nucleon). Values of A are 
chosen so that magnitudes of the resulting two-body pair¬ 
ing matrix elements are realistic. The last row of table 
U shows the average value of the pairing matrix elements 
of a specific form, mean[-Vaaaa] = - Va<i<ia/( 2 D), 

where a runs over the entire single-particle space whose 
dimension is 2D. 

The results are shown in table El We see that in gen¬ 
eral the GDM calculation reproduces the exact density 
matrices accurately, except for the very small value of 
ppp in sets. Particularly, the abrupt change of k values 


from set3 to set4 is captured very well. The off-diagonal 
parts of the density matrices (pa/i, Hap, and K-yo), 
originating from the “off-diagonal” parts of the HFB-type 
pairing Hamiltonian, are always well reproduced, ranging 
from a less-than-one-percent effect to a few. The errors 
of the GDM ground state energy relative to the exact 
ones are shown in the last column. We see that the er¬ 
rors are small and consistent with the statistical estimate 
made within a large random ensemble for the case of a 
BCS-type pairing Hamiltonian in Ref. [^ . 


V. SUMMARY 

In summary, we considered the solution of the GDM 
pairing theory with the HFB-type Hamiltonian. In the 
algorithm only the one-body density matrix is computed, 
thus the method should be much faster than the VAPNP 
-|- HFB method in which the two-body density matrix 
is needed in calculating the expectation value of a two- 
body Hamiltonian. With the assumption that the parent 
and daughter nuclei were represented by pair condensates 
with the same pair structure, the pair-transfer ampli¬ 
tudes K could be calculated in one run together with the 
occupation numbers p. In contrast, the traditional vari¬ 
ation method calculates the parent and daughter nuclei 
separately, each with a different pair structure. 

The formalism is tested in the simple two-level model 
with factorizable pairing interactions, and the semi- 
realistic model with the zero-range delta interaction. In 
both cases the GDM calculation reproduces quite well the 
exact density matrices and ground state energy within 
the physical range of parameters for the realistic spher¬ 
ical nuclear shell model. The errors are small and con¬ 
sistent with the statistical estimates made within a large 
random ensemble for the BCS-type Hamiltonian in Ref. 
( 2 ^. In the current form the theory is ready for the ap¬ 
plication to realistic nuclear systems. 
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FIG. 1: (Color online) Density matrix p in the original single¬ 
particle basis for the ensemble. In panel (a) there are 360 
points corresponding to the 360 examples in the ensemble. 
The horizontal coordinate of each point is paa of the example 
by the exact diagonalization and the vertical coordinate is 
the one by the GDM calculation. Similarly for panels (b) and 
(c), but the coordinates of the points are ppp and p^p of the 
examples, respectively. 
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FIG. 3: (Color online) Ground state energies in the ensemble. 
In the figure there are 360 points corresponding to the 360 
examples in the ensemble. The horizontal coordinate of each 
point is the pairing correlation energy impair of the example. 
The vertical coordinate Sodm is the ground state energy by 
the GDM calculation, measured from the exact ground state 
energy. 


FIG. 2: (Color online) Density matrix k in the original single¬ 
particle basis for the ensemble. In panel (a) there are 360 
points corresponding to the 360 examples in the ensemble. 
The horizontal coordinate of each point is KoS of the example 
by the exact diagonalization and the vertical coordinate is 
the one by the GDM calculation. Similarly for panels (b) and 
(c), but the coordinates of the points are and of the 
examples, respectively. 










p 



TABLE I: Particle number 2N and pairing strength A as in¬ 
put parameters for the six sets of calculations. The quan¬ 
tity b = ^h/{fiLo) in the unit of A is the length parame¬ 
ter for the harmonic oscillator wavefunctions. Also shown is 
mean[—l/Qaac], the average value of the pairing matrix ele¬ 
ments of a specific form. See text for details. 
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A/(MeV • b^) 
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40 

mean[—FQ-aaQ-J/MeV 
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FIG. 4: (Color online) Deviations of various quantities from 
the exact results as functions of the off-diagonal pairing 
strength p. The 360 examples in the ensemble are divided 
into five subgroups according to the value of p. Panel (a) 
plots the root-mean-square deviations from the exact results 
of the density matrix p by the GDM calculation within each 
subgroup, where black squares and blue circles are for paa 
and pap, respectively. Similarly panel (b) plots the root- 
mean-square deviations of the density matrix k within each 
subgroup, where black squares, blue circles, and red triangles 
are for Kqq , , and , respectively. In panel (c) the black 

squares show the average value of the pairing correlation en¬ 
ergy Epair within each subgroup, and the blue circles show 
the average value of the ground state energy by the GDM 
calculations Sgdm, measured from the exact one. 
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TABLE II: Density matrices p and k. by the exact and GDM calculations in the original single-particle basis. For clarity 
each single-particle level is represented by a Greek letter as {a,/3, 7 ,/i, = { 2 si/ 2 , 1 ^ 3 / 211 ^ 5 / 2 ) 3 si/ 2 , 2 ^ 5 / 2 }- For each set 
of calculation the results of density matrices p or k are listed in two rows, with the top row being the exact results and the 
bottom row the GDM results. The last column gives the errors of the GDM ground state energy relative to the exact one, 

EgOM = {(j>N\H\(l)N) — Eexact. 



pOLOL 

PPP 

P 77 


pvv 

Pot^l 

pju 



K'yj 




K'yi^ 

EcDivi/keV 

setl 

0.0400 

0.0536 

0.0278 

0.0313 

0.6326 

0.6261 

0.0007 

0.0008 

0.0020 

0.0016 

0.0046 

0.0056 

0.0285 

0.0227 

0.1972 

0.2283 

0.1655 

0.1760 

0.6554 

0.6524 

0.0162 

0.0171 

0.0275 

0.0282 

0.0213 

0.0224 

0.0282 

0.0227 

7.44 

set 2 

0.0889 

0.1086 

0.0342 

0.0373 

0.9460 

0.9374 

0.0009 

0.0010 

0.0013 

0.0011 

0.0078 

0.0090 

0.0233 

0.0200 

0.2910 

0.3217 

0.1819 

0.1919 

0.5837 

0.5778 

0.0178 

0.0186 

0.0265 

0.0271 

0.0243 

0.0253 

0.0136 

0.0117 

9.95 

sets 

0.9766 

0.9423 

0.0420 

0.0312 

0.9776 

0.9968 

0.0030 

0.0029 

0.0012 

0.0007 

0.0479 

0.0454 

0.0216 

0.0183 

0.9418 

0.9041 

0.1482 

0.1762 

0.1932 

0.2148 

0.0276 

0.0283 

0.0163 

0.0178 

0.0449 

0.0423 

0.0039 

0.0036 

129 

set 4 

0.9867 

0.9934 

0.5254 

0.5085 

0.9851 

0.9947 

0.0030 

0.0025 

0.0014 

0.0010 

0.0478 

0.0478 

0.0204 

0.0189 

0.1419 

0.1277 

0.7048 

0.7000 

0.1455 

0.1367 

0.0141 

0.0138 

0.0249 

0.0248 

0.0062 

0.0055 

0.0024 

0.0021 

62.9 

sets 

0.9952 

0.9935 

0.0078 

0.0029 

0.9959 

0.9999 

0.0007 

0.0007 

0.0002 

0.0001 

0.0231 

0.0228 

0.0096 

0.0089 

0.9930 

0.9914 

0.0519 

0.0540 

0.0576 

0.0586 

0.0124 

0.0127 

0.0046 

0.0046 

0.0227 

0.0225 

0.0005 

0.0005 

42.1 

set 6 

0.9032 

0.8407 

0.2134 

0.2428 

0.8775 

0.8795 

0.0139 

0.0131 

0.0079 

0.0074 

0.0990 

0.0931 

0.0541 

0.0471 

0.6378 

0.5500 

0.4148 

0.4623 

0.4880 

0.5086 

0.0588 

0.0587 

0.0652 

0.0708 

0.0641 

0.0553 

0.0255 

0.0236 

76.7 




